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Abstract 

We describe a possible search for muonic radioactivity from lead nuclei 
using the base elements ( "bricks" composed by lead and nuclear emulsion 
sheets) of the long-baseline OPERA neutrino experiment. We present 
the results of a Monte Carlo simulation concerning the expected event 
topologies and estimates of the background events. Using few bricks, we 
could reach a good sensitivity level. 



1 Introduction 

In the late 80's, some theoretical work was dedicated to investigate possible 
exotic types of nuclear radioactivity, consisting in the emission of light particles 
such as pions or muons from heavy nuclei [1, 2]. Although the theoretical esti- 
mates indicated a larger branching ratio for spontaneous (or neutron-induced) 
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muon emission than for pions 1 , the experimental searches were mostly dedi- 
cated to the pionic radioactivity [3, 4] , since the detection conditions were more 
favourable with the used experimental set-ups. None of the experiments us- 
ing heavy radioactive nuclei could prove the existence of the pionic or muonic 
radioactivities, but some of them may have produced indications in support 
of these hypotheses [5] . Some experimental indications for pionic radioactivity 
come from the interpretation of radioactive super giant halos in crystals [6]; 
such a technique, however, cannot give precise information concerning the na- 
ture of the source or allow quantitative estimates of branching ratios, lifetimes, 
etc. A short review of the theoretical problems and of the experimental results 
obtained, mostly for the pionic radioactivity, may be found in rcf. [7] . 
Muons or pions could be emitted by nuclei through the decays [2] : 

(A,Z) ^ fJ , ± + ^(V^) + (A 1 ,Z 1 ) + ... + (A n ,Z n ), (1) 

{A, Z) -> tt^tt ) + (A u Zi) + ... + (A n ,Z n ), (2) 

where, for reasons of energy and momentum conservation, the number of frag- 
ments n is > 2. The nuclei (A\,Z\), (A 2 , Z 2 ), {A n , Z n ) would yield a se- 
quence of (5~ decays leading finally to stable nuclei with a balanced neutron- 
proton ratio. 

In rcf. [2] some nuclear charge thresholds for different possible spontaneous 
particle emission were listed: 

(i) fi^ (prompt muon) for Z > 72 

(ii) tt^~ (prompt pion — > delayed muon) for Z > 76 
(Hi) 2^i ± (prompt muon pairs) for Z > 91 

(iv) 27r ± (prompt pion pairs — > delayed muon pairs) for Z > 100 

Natural lead is mainly composed by three nuclides: 206 Pb (24.1%), 207 Pb 
(22.1%) and 208 Pb (52.4%). They are stable nuclides, but the channels (i) and 
(ii) are energetically allowed. 

The spontaneous or neutron induced fission of Pb has never been observed. 
Nevertheless, the decay into a lepton-antilepton pair (muon and neutrino) or 
into a quark-antiquark pair (forming a pion) could enhance the probability that 
the remaining "hyper-cold" nuclear state be unstable [1]. 

A search for a Pb muonic decay can be made as a byproduct of the OPERA 
experiment [8]. Uranium or Thorium could be used in the future; such possible 
experiments would benefit from the results of an initial search using Lead. 

In the hypothesis of the decays (1) and (2), the fission fragments would 
remain nearly at rest; most of the available energy would be used to produce 
the /j, (or 7r) and the kinetic energies of [i and (or it). 

The total kinetic energy in a decay 

Pb -> / i ± + V fl (pn) + (A l ,Z l ) + (A 2 ,Z 2 ), (3) 

^^This observation is not necessarly valid in the case of Pb, as the predicted branching ratios 
for the pionic and muonic radioactivities are defined in [1, 2] relatively to the spontaneous 
fission of the parent nuclei. 



2 



assuming close values of A\ and Ai (symmetric fission), is about 30 MeV for 
negative muons and about 20 MeV for positive muons. Considering that the 
associated muon neutrino takes away a sizable fraction of this energy, the spec- 
trum of emitted muons could be like in a (3 decay, with an average around 10-^15 
MeV. 

The decay into two fragments would not lead to the best energetic situation 
since the two resulting nuclei would on the average have large atomic mass and 
too large neutron numbers. Decays into more than two nuclear fragments would 
lead to more energetic muons; at the moment there are no estimates available 
for such decays. 

There is only one experiment that published upper limits for exotic decays 
of Pb [9]. The experiment used 123 g of natural lead, obtaining a 90% C.L. 
upper limit on the 7r° — ► 77 counting rate of 3.3 • 10~ 27 s _1 , for a total counting 
time of 109 hours. The same authors searched also for pionic emission from 
Uranium, and reported a 90% C.L. upper limit for T^/Tsf ~ 3.1 • 10" 4 , where 
Tsf is the width of the spontaneous decay fission. Using 235 U and 252 Cf they 
obtained upper limits of 1.4 • 10 -4 and 3.3 • 10~ 10 , respectively. A search for 
spontaneous emission of both muons and pions from 252 Cf yielded upper limits 
for the ratio T^/Tsf in the range 10~ 6 10~ 8 [3]. 

There are no experimental limits for muonic radioactivity. 

In the following we describe a possible search for muonic radioactivity of 
lead nuclei and we discuss the experimental set-up, the results of a Monte Carlo 
(MC) simulation giving estimates of the geometric efficiencies, a description of 
possible backgrounds and the reachable limits. 

2 Experimental set-up 

We propose to perform an experimental search for muonic radioactivity from 
lead nuclei in the low background conditions offered by the Gran Sasso under- 
ground Laboratory (LNGS). The very low cosmic muon flux and the low natural 
radioactivity of the rock in the experimental halls of the LNGS provide unique 
conditions, allowing a potential discovery, or at least to establish a very good 
upper limit for this exotic decay process. A detailed description of the different 
background sources is given in Sec. 4. 

This search could be a by-product of the OPERA experiment [8] presently 
under construction at LNGS. OPERA is a hybrid long-baseline neutrino ex- 
periment, aimed to the direct observation of v T appearance from < — > v T 
oscillations in the CERN-Gran Sasso v„ beam. The detector is made of 2 super- 
modules, each with a massive lead/nuclear emulsion target, electronic detectors 
and a magnetic spectrometer. Nuclear emulsions are used as high resolution 
(0.06 ^m) tracking devices for the direct detection of the decay of the r leptons 
produced in the charged current v T interactions with the target. 

The OPERA base element ("brick") is composed of 56 lead sheets (1 mm 
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Figure 1: Illustration of the sequence lead/nuclear emulsion sheets which char- 
acterise an OPERA brick. 



thick) interleaved with 2 nuclear emulsion sheets (43 /<m thick) on both sides of 
a 200 /im thick plastic base. This sequence is illustrated in Fig. 1 . The area of 
each sheet is 10.3 x 12.8 cm 2 . 

The OPERA lead bricks (each containing a mass of 8.23 kg of Pb) could 
allow an experimental search for muon emission from lead, with exposures of 
several months. Their analyses with the fast automatic optical microscopes [10] 
would establish the local background contributions and validate the analysis 
procedures. 

A preliminary test using a "small" brick composed of 22 lead sheets and 
21 recently refreshed nuclear emulsion sheets was started at the end of March 
2005. As the background rejection/reduction (see Sect. 4) is a crucial point 
for this search, the detectors are surrounded on all sides by a shield, making a 
closed box structure similar to those used in experiments for dark matter and 
for neutrinoless double beta decay searches. The shield is composed of an inner 
layer 5 cm thick of very pure copper followed by 15 cm of very low activity 
lead 2 . The third and last layer of the shield is a 3 cm thick polyethylene, in 
order to absorb neutrons. A scheme of the installed shield, with the dimensions 
for containing 5 OPERA bricks, is presented in Fig. 2. The set-up is located 

2 These thicknesses have been found to be adequate with a Monte Carlo simulation con- 
sidering the effects of 0.5 4- 2.6 MeV photons, which could produce electrons mimicking the 
searched events. 
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76 cm 



Figure 2: Sketch of the shield used for background rejection/reduction. It 
is composed by layers of very pure copper and lead, plus an outside layer of 
polyethylene ("transparent" box). The set-up is located in the OPERA emul- 
sion storage room, equipped with a ventilation system. 



in the emulsion storage room, in hall B of the Gran Sasso Laboratory. The 
radon reduction is obtained with a ventilation system due to fresh air forced 
circulation. 

The mean range of 15 MeV muons in lead is about 1.5 mm, so most of the 
possibly produced muons could be detected in the two layers of emulsion close to 
the Pb sheet "source" ; depending on the angle of emission and on the location 
of the decay point, some of them could also be detected by the next emulsion 
sheets. 

It is clear that thinner lead sheets would allow the emitted muons to loose a 
lower fraction of their energy to escape from the lead and to cross more emulsion 
layers. They would also yield a higher detection efficiency. 

3 Monte Carlo simulation 

A MC simulation program was implemented to estimate the occurrence of dif- 
ferent event topologies. The simulation is based on the GEANT3 [11] package 
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Figure 3: Distribution of the number of emulsion layers crossed by the simulated 
fi~ (left) and by the electron resulting from its decay (right), assuming an initial 
muon energy of 15 MeV and isotropically distributed emission angles. 



applied to the OPERA lead/emulsion set-up. At this point it does not include 
the complete detector response and does not estimate the real event reconstruc- 
tion efficiencies. 

The simulation reproduces one complete OPERA brick, where muons of 
different energies (see the first column of Table 1) originate in random positions 
in the lead sheets. The initial muon directions are isotropically generated. We 
assumed different definitions for a candidate event, requiring that the muon 
crosses at least: (i) one single emulsion layer, (ii) two emulsion layers (near the 
same base), and (iii) three emulsion layers (and thus also a lead plate, which 
would allow a better identification using the measured energy loss). We also 
requested the detection of the decay positron or electron in at least (iv) two, 
(v) three, or (vi) five emulsion layers, together with the muon detection. 

Fig. 3 shows the distribution of the number of emulsion layers crossed by 
a sample of simulated \i~ (left) and by the electrons resulting from their decay 
(right), assuming an initial muon energy of 15 MeV and isotropically distributed 
emission angles. 

As there is no theoretical prediction concerning the muon energy spectrum, 
an overall fraction of events with a defined topology cannot be estimated. The 
percentages of events listed in Table 1 were computed for samples of MC events 
with fixed energies. The detection threshold for the [i ± , obtained requiring at 
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(a) (b) 

Figure 4: Simulated events of spontaneous [i + (a) and ji~ (b) emission from lead 
inside an OPERA brick, assuming an initial muon kinetic energy of 15 MeV. In 
(a) the /i + starts from the second lead sheet and stops in the third one, where 
it decays into e + v e V^. The fi + is seen in 2 emulsion sheets, the e + in 4. In (b) 
the \x~ starts from the fourth lead sheet and stops in the second one, where it 
decays into e^v e v^. The \i~ is seen in 4 emulsion sheets, the e" in > 6. The 
change of direction for the /z and in the interfaces lead-emulsion are due 
to the MC approximation which, for each material, adds the multiple Coulomb 
scattering effects and gives a global scattering angle of the outgoing direction 
relative to the incoming one. 



least 2 nuclear films crossed by the ^ ± , is about 3 MeV. For each topology, the 
values quoted in Table 1 were obtained as averages over all possible muon emis- 
sion directions. These estimates may be considered as geometrical efficiencies 

In Fig. 4 are shown two simulated muon emission events from the lead 
inside an OPERA brick, assuming an initial kinetic energy of 15 MeV, and the 
production of a //+ in Fig. 4a and a \i~ in Fig. 4b. The dashed lines (longer 
lines) represent the muon trajectories, the solid lines are positrons or electrons, 
the dash-dotted lines are neutrinos and the dotted lines are photons. In this two 
favourable cases the muon and the electron/positron would be easily identified. 

3.1 Estimates of global detection efficiencies 

The geometrical e g have to be multiplied by the and e^- reconstruction effi- 
ciencies to obtain the total efficiency. We made a first estimate of ttot for muons 
based on the measurement procedures and the algorithms used in OPERA [10]. 
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Table 1: Percentage of events as a function of the minimum number of emulsion 
films crossed by both \i (N^) and e (N e ), for different [i initial energies (col. 1). 
The column in bold refers to the conditions requested by the OPERA tracking 
system to define a track (N^ — N e = 4). 



The relation between the energy of charged particles and the density of 
grains produced in emulsions along their trajectories [12] shows that muons 
with the predicted energies would give recognisable black tracks without relevant 
scattering in the emulsion films. For the decay electrons/positrons the tracks 
would be similar, but with a lower density of grains per unit path-length. 

With the OPERA tracking procedure [13], the mean detection efficiency for 
each "microtrack" in one emulsion film is ~ 95%; it is a value averaged on the 
angular range between tracks normally incident on the plate surface and tracks 
with incident angles ~ 0.8 rad with respect to the normal. The instrumental 
limit of 0.8 rad on the incident direction of each microtrack introduces an event 
selection of ~ 30%. 

The "base track" (obtained from 2 microtracks separated by the plastic 
base) reconstruction efficiency of ~ 90% (95% x 95%) is comprehensive of the 
uncertainties due to the emulsion shrinkage and possible surface distortions. 
The top-bottom linking efficiency ranges from the ~ 50% for 5 MeV particles 
to ~ 99% for particles with energies > 20 MeV. It is mainly due to the multiple 
scattering in the plastic base. The base track linking efficiency is ~ 6% for the 
whole range of initial muon energies. 

In order to reduce the background, an algorithm defines a "track" requiring 
at least 4 emulsion films (it means at least 1 lead sheet and 2 plastic bases) 
crossed by a particle. The detection efficiency e to t for a muon crossing at least 4 
emulsion films is the product of the percentage of events with these topological 
requirements e g (col. 6 of Table 1), the 30% given by the angular limit, 4 times 
the 95% efficiency for the microtracks, 2 times the top-bottom linking efficiency 
and the linking efficiency between 2 base tracks. As shown in Table 1, the 
threshold muon energy to have at least 4 emulsion films crossed by the muon is 
10 MeV. 
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For the electrons, the probability of crossing at least 4 emulsion films e e g is 
about 68% for the whole range of initial muon energies. This percentage was 
obtained from the simulated events requiring a minimum number of 4 emulsion 
films crossed by the electron and no constraints on the muon trajectory. The 
convolution with the topological requirements of the angular directions within 
0.8 rad (30% selection), the 95% efficiency reconstruction of the 4 microtracks, 
the top-bottom linking efficiency and the linking efficiency between base tracks 
yields the total electron detection efficiency. 
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Table 2: Geometric efficiencies (e g ) and total muon detection efficiencies 
(e^ f ERA ) as a function of the initial muon energy computed on the basis of 
the OPERA tracking procedure for — 2 (col. 2-3) and for = 4 (col. 
4-5). The values t* tot (col. 6) are obtained relaxing the angular and position 
tolerances in the base track reconstruction (see text). All values are in %. 

We intend to optimise the tracking procedure for this search using new scan- 
ning conditions, ft may be possible to require only 2 emulsion films crossed by 
a muon to reconstruct its track. This optimisation should be possible, because 
of the easily recognisable black tracks produced by slow muons. Once a muon 
is found, we should be able to check visually the presence of the decay electron. 

Table 2 shows the % efficiencies e g and e^ Q f ERA computed requiring at least 
2 emulsion films crossed by the muon (col. 2-3) and 4 emulsion films crossed 
by the muon (col. 4-5). The values refer to the present OPERA tracking 
[13]. The linking efficiencies may be improved relaxing the angular and position 
tolerances between consecutive base tracks. This should be possible because of 
the low background conditions. Requesting an angular tolerance of 100 mrad 
and a position tolerance of 100 fim, we obtained the values given in col. 6 of 
Table 2. 

In OPERA the different emulsion sheets of each brick are aligned mechan- 
ically and then exposed to cosmic ray muons to have a relative alignment of 
~ 1 (Um [14]. In the present case, due to the low background conditions, we 
expect that the mechanical alignment precision should be adequate for the pur- 
pose; since the emulsion layers are placed horizontally, a check of the alignment 
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could come from the few high energy cosmic ray muons which should cross the 
emulsions [15]. 

As the expected halflifes t\/2 are much larger than any reasonable exposure 
time T, the expected sensitivities are estimated from 

ln2 m 

N~ = l— T£ tot (4) 

where SN — 2.3 is the number of events corresponding to a 90% C.L. limit 
assuming no candidates, N is the initial number of nuclei, and e tot is the ex- 
perimental efficiency 

Assuming the use of one OPERA brick for one year exposure and a global 
detection efficiency of ~ 10%, we could reach a sensitivity of about 7 • 10 23 yr 
(90% C.L.). 



4 Background estimates 

The background originates from several different sources. 

1) The environmental radon background is reduced by the ventilation and 
may be monitored using nuclear track detectors, such as CR39 or Makrofol, 
that are insensitive to muons. 

2) The background produced by the ambient neutron flux (~ 0.42 • 10~ 6 
neutrons/s/cm 2 , with energy greater than 1 MeV [16]) should be studied; its 
effect can be reduced by a proper shielding and appropriate "event definitions" . 

3) The background produced by the radioactive 210 Pb isotopes present in 
the lead. The a and (3 radioactivity of the lead used in OPERA [8] is < 0.02 
Bq/g. A single lead plate has an activity < 3 Bq, and only a fraction of the 
particles emitted from the outer surfaces of the material escapes from the lead 
sheets. The a background could be eliminated by thin plastic foils inserted 
between the lead plates and the emulsion sheets. 

4) The background produced by radioactive nuclides present in the emulsion, 
for example 212 Po nuclei. The emitted a particles have an energy of 8.785 MeV 
and a range in emulsion of 74 fxm. This background should be reducible by 
dE/dx measurements on the tracks and from range considerations. 

5) The background due to cosmic rays. 

(i) Cosmic ray (CR) muons reaching from above the LNGS underground 
labs (3700 hg/cm 2 average rock cover) are ~ 1 muon/h/m 2 [15]; the number of 
CR muons expected to cross one brick during one year is about 90, sufficient 
to improve alignments of the different emulsion layers. They have an average 
energy of 240 GeV [15]. Tracks produced by CR muons in emulsion may be 
easily identified on the basis of their high energy and relatively small energy 
loss. Low energy CR muons can be removed requesting that the tracks start in 
a lead plate. 
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(ii) Pions produced by muons interacting in the lead plates. Slow pions 
decay into muons and simulate a spontaneous muon emission. This process was 
investigated by MC techniques in MACRO, since it was a background for the 
measurement of the upgoing muon flux [17]; the probability to produce a pion 
by a CR muon through photonuclear interactions in a layer of 150 g/cm 2 was 
estimated to be < 2 • 10~ 5 , and the fraction of low energy backward scattered 
pions is < 2 • 10~ 4 . Since in an OPERA brick the yearly rate of charged pions 
produced inside the plates is < 6 • 10~ 4 , this background may be neglected. 

(Hi) The flux of muons resulting from atmospheric interactions around 
the detector is known from measurements done by MACRO [18]; they yield a 
negligible background. 

(iv) In the Gran Sasso underground laboratories the CR neutron flux is 
about ~ 10 neutrons per square centimetre in one year [16, 19]. Only a fraction 
of these neutrons interact in the detector and the events produced by them 
should be recognisable in emulsions. 

6) The future CERN-Gran Sasso neutrino beam may yield a background 
contribution; the energies of muons induced by the v beam will be much higher 
than the energies of emitted muons from lead. Moreover, the beam will be 
continuously monitored. Also this background should be easily recognised and 
removed. The beam will not be present before the middle of 2006. 

5 Conclusions and perspectives 

We described a test search for spontaneous emission of muons from Pb nuclei, 
using some OPERA lead/emulsion bricks. The proposed test is suggested by 
some theoretical works performed in the last years; it would be complementary 
to other experiments looking for such exotic radioactivities from heavier nuclei 
and could reach a good sensitivity level. For a visual fine-grained technique (the 
emulsion technique) the identification is easier and the background lower. The 
test will yield the first information relative to the background in the OPERA 
bricks. 

We are in the process of making a complete simulation of the detector in- 
cluding its response and the track reconstruction efficiencies. As already stated, 
thinner lead sheets would improve the scanning and the track reconstruction 
efficiencies. 

We have shown that stringent limits for spontaneous muon radioactivity may 
be reached. We would obtain ii/2 > 7 • 10 23 years. 

In the future, one could repeat the experiment using elements heavier than 
lead. 
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